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ABSTRACT 


Consider the two-factor crossed components of variance model given 
by 

Y iJ = y + A i + T J + E iJ for 1 = 1 ’ "- 1 * J = 1 . J * 

The random variables A^, Tj, are normal and independent with means zero 

and variances a 2 , a^, a 2 . Approximate confidence intervals are presented and 
evaluated for the ratios p, * a 2 /(o 2 + c 2 + o 2 ) and p_ = ar 2 /(a 2 + a 2 + a 2 ). 
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1. Introduction . 

Consider the twu-factor crossed Components-of-Variance model 
given by 

hr W £ « 1=1.2....1; J-1 >2.--J • 

The «‘*'e independent unobservable random variables and 

A i ~ *M(0 ,o a 2 ) Tj ~ N(0.o t 2 } - N(0,a e Z ) 

H is an unknown parameter, and the are observable random 

variables. 

The analysis of variance table for this model is given in Table l. 
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— ? 2 2 

The random variables Y, S.j Sg » are complete sufficient 

2 ? 

statistics for this problem and /e^~ x (n^for i=1,2,3 and 
are mutually independent. 

An important problem in applied statistics is to obtain confidence 
Intervals on various functions of the variance components o^ 2 , o.j. 2 , and 

For some well-known functions exact-size confidence intervals a*e 
available, but in mar.y situations exact-size confidence intervals have 
not as yet been derived. 

For the following important and useful functions exact-size 
confidence intervals are known (exact-size means exact confidence 
coefficient): o e 2 , ° 2 /q£ * o_ Z /o T 2 , 0 L 2 /(° C 2 4 ° A 2 ). ° e 2 /(° e 2 4 °j 2 )» 
o a 7(o £ Z + o^ 2 ), Oj 2 /(u e 2 + Oj 2 ). For a discussion of these see Greybiii (1976). 

O 

Very good approximate confidence intervals are known for o A 
and Oj 2 

For approximate confidence intervals on linear combinations of 
2 2 2 

°A ’ °T * °e see We ^ c * 1 (1956). 

In applied work it is often important to determine the ratios of 
2 2 2 

the Individual variances , and o £ to the total variation 

2 2 2 

( a^ + Oj + o ), i.e., to determine 

« 2 y / 2 . 2. 2i 

P A °A /(o A + °T + ° c } 

. Z,l 2 . 2 . 2 . 

Pj ■ Oj /(o A + Oj + 0 £ ) 

2,r 2 . 2 . 2, 

p * o /(o A + o T +o ) 
c c “ • e 
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There Is no method available for setting exact 1-a confidence 
Intervals on p^ , py , or p £ . Large sample approximations have 
been proposed by Osborne and Paterson (1952). 

The purpose of this article is to obtain good approximate 1-a 
upper, lower, and two-sided confidence intervals for p^ and Py 
and evaluate how close the approximate confidence coefficients arc to the 
specified 1-a . 

The quantities p^ and py are useful In many applications. 

Graybill (1976) gives an example of a two-factor crossed components of 
variance model without interaction. An investigator wants to examine 
the effects of drivers and automobiles on gas mileage. One type of 
automobile is used and 20 are selected at random from one month's 
production. Ten people are selected at random from a large city and 
are asked to drive each car for two weeks and record the number of miles 
per gallon of gas used for each car for the two-week period. The total 
variation In the number of miles per gallon of gas used for this type of 
automobile is + °f + a \' w * iere °j| tlie variation due to the cars, 
o T is the variation due to the drivers, and o ? is the variation due to 
all other noncontrolled factors. The proportion of the variation in the 
number of miles per gallon due to the cars is given by , and the 
proportion of the variation in the number of miles per gallon duo to the 


drivers is given by py. 






2• Approximate C o nfidence Intervals . 

In this section approximate confidence intervals on p A and 
Pj will be derived. 

First consider a 1-a upper confidence interal on p A 
denoted by < o A where is the lower confidence point. 

It is simpler to find a 1-a upper confidence Interval for 

Jb fl 2 

e. » 


Ko t z + o c 2 ) 


e l ‘ 9 3 
e 2 + "l e 3 


and convert this to a confidence interval for p a by using the fact 
that / \ 

p a a yrr-fe-J (2 - 1} 


A function g L (Y) of the observable random vector must be 
determined such that 

P[9 L (Y) £ 0 A ] 

Is close to a specified 1-a for all values of the unknown parameters. 
2 o o 

Since Y, , S 2 . S 3 is a set of complete sufficient statistics 
for this problem, the lower confidence point should be a function of 
these random variables. So a function h L (Y,S^ ,Sg ,S 3 ) must be 
determined such that 


P[h L (V.s 1 2 .s 2 2 ,s 3 2 ) < e A ] 


Is close to a speclfieo 1-a . 

If any constant k is added to each observation Y^ the 

2 2 2 

parameter o A is unchanged so it is required that h^(Y+k,S^ ) 
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remains unchanged. The constant k = - Y is chosen and thus 

h L (T,S 1 2 ,S 2 2 ,S 3 2 ) = h L (0,S 1 2 ,S 2 2 ,S 3 2 ) = q L (S 1 2 .S 2 2 ,S 3 2 ), which is a 
2 2 2 

function of ,S 2 , and S 3 only. 

If each observation Y^j is multiplied by a non-zero constant 

c , the parameter 6^ is unchanged so it is required that 
2 2 2 

^(S] ,S 2 »S 3 ) remains unchanged also. The constant chosen is 

1/S 3 2 so / \ / > 

Is 2 s 2 \ U 2 s 2 

V,S 3 2 ) - %[ ^ ^ • l] - L ^2 * S 


2 2 2 2 

which is a function of /S 3 = F^ and S 2 /S 3 = F 2 only. 


The problem is reduced to determining a function L(Fj,F 2 ) such 


that 


P[L(F r F 2 ) < e A ] 

will be close to a specified confidence coefficient 1-a . 

To determine the function L(F^,F 2 ) the following five conditions 
are imposed. 

1) The confidence interval must be exact when e 3 = 0. 

2 

When e 3 = 0 it follows that S 3 = 0 with probability one, 

Fj~^ ■ 0 with probability one, and = 0j/e 2 ' an exact 
lower confidence limit on 0^ is 



F 



2 

2 


and L(F^,F 2 ) must be equal to 








7 


a:n-| ,n 2 2 


In order to satisfy this condition. 

2) When o ^ -*• it is required that the confidence coefficient 
2 2 2 

be exact. When ■** and Oj ,o £ are fixed it follows that 

6^ -*® with © 2* e 3 fixed, and dominates e A . So the confidence 
Interval on should behave like an exact 1-a confidence interval 


on 


6j . So L(Fj,F 2 ) should be equal to S^F^" 1 


3) The confidence interval is required to be exact when 
J -*» and I is fixed. 

2 

When J -*» it follows that n 2 •*»> n 3 -*«> and hence S 2 -*■ e 2 
In probability and S 3 2 -*■ e 3 in probability; also, 

(S 2 2 + n^S, 2 ) - (e 2 + 0,63) 

In probability. 

Begin with an exact (1-a) confidence interval for e 1 given 
by 

s , 2 f “„ < e, 

and subtract S 3 2 from the left side and the "equivalent" value e 3 
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from the right side to obtain 

S 2 F - S 2 <9,-0, 

1 3 — 1 3 

Then divide the left and right sides, respectively, by 
2 2 

(Sg + n j S 3 ) and the "equivalent" value (e 2 + n^) to obtain 
s 2 F - 1 - s 2 

1 a:n,,« 5 3 e, - 9. 


’1 °3 


c 2 , p 2 6, + n,9, 

S 2 + n-jS 3 2 1 3 

Hence, when J-*» the lower confidence point should behave like 

L ( f r F 2 ) ‘-r-frp 

4) When I-*° and J is fixed the confidence interval is 
required to be exact. 

For I-*- and J fixed it follows that n^-*® , n 3 -x» and 

2 2 
S-j -*■ 9^ in probability, and -» 9^ in probability; also, 

2 2 

S| - e i " e 3 in probability. 

Begin with a 1-a exact confidence interval for e 2 given by 

09 < S 2 F . „ 

2—2 a:»,n 2 

2 

Add n-|S 3 to the right side and the "equivalent" value n^e^ to the 
left side to obtain 

9 n + n,9, < S, Z F . + n. S» 2 

Z 13—2 a:»,n 2 1 3 

Take the inverse of each side, then multiply the right and left sides 
respectively, by S-| - and-its "equivalent" value o-j - 
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to obtain 


S 2 - S 2 


s o F . „ + n,S * 

2 a:«,n2 1 3 


9 1 - 9 3 
- e 2 + n^ 


Hence, when I-*~,L(F^ .Fg) is required to behave like 


L(F r F z ) = f 2 f J — + n ' 

2 a:«,n 0 1 


5) If an a-level test of H Q :e^ = g 3 vs. > e 3 is 

accepted the confidence interval should include "zero" and L(F^,F 2 ) 
should be increasing in for fixed values of F 2 . 

H is accepted iff F, < F „ , so the confidence limit 

o i — a.n>j ,n 3 

should be 

l (F,,F 2 ).0 when F, < F^ ^, 

L(F r F 2 ) monotonic increasing in F-j for fixed values of F 2 

when F, > F _ 

I a:ni»n 3 . 

The form of the function that could be used for L(F^,F 2 ) is 

a„ + ^iFt 

L(F 1 ,I V = b + b,F ? 

0 1 2 

but a more general function will be examined, i.e., 


L(F r F 2 ) 


.V*lV a 2 f l 


b Q + b,F 2 + b 2 F 2 ‘ 


where a., and b^ which are functions of n^.n^n^.r* are determined 
by requiring L(F 1 ,F 2 ) to satisfy conditions 1), 2), 3), 4), 5) 
Substituting condition 1) gives 


± = f r 1 


(2.2) 
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Substituting condition 2) gives 


and hence 


a, = F 

1 ain-j," 


^1 ^a:n,,»^a:n, ,n. 


i-| 1— u. ■ i-| >■ >2 


(2.3) 


(2.4) 


-1 


From condition 3) Lim a„ = -1 , Lim b = n, , Lim a, = F „ 
J-*« 0 J-*» 0 1 J— 1 a: V 


Lim b. =1, Lim a 2 = 0 , Lim b ? = 0 

J-Ko J-KO J-X» 


which are satisfied for those values of a-|, b^ from equations 
(2.3) and (2.4). 

From condition 4) Lim a = - I, Lim b = n, , Lim a, = 1 » 

I-XX, 0 I-X» 0 1 I-KX. 1 

Lim b, = F „ , Lim a 0 = 0 , Lim b, = 0 . 

1 a:», n 9 , 2 r, 2 

c 1'Ko 

Hence the limits of and b^ are consistent for condtions 3) 
and 4). 


From condition 5) L(F 1> F„) =0 if F, < F _ which implies 

'12 1 — a: n-j,n^ 


that 


a rt + a i F i + a 5 F i" = 0 if F, < F . 

o li 2 1 I — a:n.j ,n^ 


Now substitute F, = F _ into 

I a .r»i » n 3 


■i 


a„ + a,F „ + a,F „ = 0 , 

0 1 C(• 0| 2 a.n-| , 

let a-j be the value given in equation (2.3) and solve for a 2 to 

obtain 


a 2 ^ a o ^a:n^ ,<»^a:n^ .n^^ain^ ,n^ 
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Use the simplest values for a Q and b Q , namely, a Q = -1 

and b ® n, , and obtain 
o.l 

- 1 + Ml,,- F 1 + (1 - F all, - F a:n 1in > F a : ni n, F 1 

L( F F? )- } - - T -^ 

n, + F "I F _ F, + b„F 


-1 


1 am.,,® a:n^ ,n^ 2 2 2 


if F i > F n n 
1 a H|,nj 


(2.5) 


L(F r F 2> ' 0 ff F li F . : n,,n ; 


To determine the value of b^ we Impose condition 0). 

2 2 

Condition 6;, When H Q :oy =0 vs. H^Cy >0 is accepted 

the confidence interval is required to be exact. 

2 2 

To examine condition 7) we get H Q :oy =0 vs. H^:oy >0 <“=> 
H 0 :9 2 = 8 3 vs. H,:e 2 > e 3 <-. H o :0 a - (8,-e 3 )/Io 3 vs. 

H l :e A < (°r«3 )/l9 3 • 

If H is true then 8. = Jo.^/Io ^ 

O A A e I 

So if H q is accepted at a level we want L(F^,F 2 ) to be an 

exact 1-a lower confidence point for . 

H is accepted iff F„ < F „ _ . An exact 1-a lower 
o 2 — a:n2in2 

confidence point for - 1/1 is 


s 2 F - 1 

1 a;n l ,n 3 _ l 9 1 ~ 9 3 
IS, 2 1 “ 1 °3 


ft-) 


So when F, < F . 


2 — a:n 2’ n 3 we want 


s 2 f - 1 

1 a:n 1 ,n~ 1 

L(F r F 2 > ‘ 77z 


IS. 


I * 


This is not consistent with the firii. live coik. : t io.' . so we only require 
the denominator of (2.5) to be I when condition <<) is satisfied. 
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So set 


and then 


1 ci:n -j a:n^ .n^ oin^in^ 2 oiin^tn^ 


-1 


^2 ^ ^a:n^ ,“^a:n^ ,02^0:ngtn^^a:in^ 


and this satisfies Lim = 0 

J-*vO 

The final result is 


-1 + F F, + ( 1 -F “If. )F f: 

a:^,®® 1 ' a:n-|ain^.n^' arn^n.^ 1 


Uf^fJ •-n- 1 -H-^^-rr 

* c n+F* F f +M-F F F \ r F* 

1 a:n| j®’ a:ni ,n^ 2 a:n^,“ an-j .n^ a:n2»n2 ain^.n^ 2 


If 


j S l 2/S 3 2 i f a: ni ,„ 3 

and 

S? 2 /So 2 > F 
l o ci • rig $ 


-1 + F :I F, + (l-F F _ )F n f: 1 
L (F r F 2 )= -3- 1 -3—1-3_J- (2.6) 


if 


s2 /S3 2 > F 

1 3 a.rij 

and 
2 2 

So /So < F„. n n 

2 3 a.r^iHs 


LtF,.F 2 ) = 0 if V/i/if.:,,,., 
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Instead of using L(F-j .F^) in equation (2.6), the simpler function 
in equation (2.7) will be used; which is obtained from equation (2.5) 
by setting b 2 = 0 . Note that Lim b 2 in equation (2.6) satisfies 

1-Kd 

the requirement in condition 5). So the lower confidence point for is 


L A^ F 1’ F 2^ 


-i + r tI f, + (i-F _1 f )f f ; 1 

a.n-j,* 1 a:n^ a:n i* n q a:n i ,n 3 ^ 


n l F a:n^,« F a:n^,n 2 F 2 


if S 2 /S 2 > F 

1 3 ex • n<j)rij 

2 2 


(2.7) 


L A (F r F 2 ) = 0 if V^F^^ 

Note that the confidence interval in equation (2.7) covers the 
maximum likelihood estimator of 0^ , say o A . To demonstrate this, 
the M.L. estimator of , denoted by 9^ is 


e A = (Ft- 1) /(F 2+ n !) if F ] > 1 


e A * 0 


if F. < 1 


and the confidence interval is 


-1 + F F, + (1-F F )F F “ 1 

a.*n^,» l a:n^ ,<« ain-j.n^ a:n^,n 3 1 


n, + F ^T F -n „ F 9 
I a. n-|,“ a:n-j ,2 


— 0 A <• 


0 < 6 . <~ 
- A 


if F i > F n 
1 oi.n^ .n^ 

1f F 1 < F . n n 
I — a: n^ .n^ 








u 


-1 + F. 


t:n-j»— F 1 + /° F a:n 1 ,n 3 ^a:n 1 ,n 3 F l < F 1 ‘ 1 


F ' < 1 , 

a:nj,« 


(1-F F . ) < o. 

o• ni f°° a.n 1 j 


n i + F "][ F . F, > F, + n, 
I a:n^,» a:n^,ng 2 2 1 


because 


a:rij ,«o a:n^ .n^ 

So L A^ F l ,f V < e A < “ anc< e A * s in the l" a confidence Interval. 

Next consider a lower l-o confidence interval for e A . A lower 
1-a confidence interval for o A> given by 0 <_ o A 5 U A^ r i* F 2 ^ can be 
readily obtained by replacing a by 1-a in equation (2.7). The 
upper 1 -a confidence point of the lower confidence interval is 


given by 


u A ( F r F 2 ) 


l-a:n.|,oo F l + ^ F l-a:n^ ,® F l-a:n, ,n_^ F l-a:n, ,n F 1 


l U.ll<| I "a • BI «| j 

n, + F,F, F_ 

I 1-a:ni,® l-ain^,n^ 2 


U A (F^» F 2 ) = 0 


If F i > F i •„ n 

I l-a.n^ * 

If F 1 1 F 1 n 

1 I-a:n^ ,n^ 


The procedure for determining this lower confidence point is 
similar to that described for the upper confidence point. 

From the formulas for upper and lower confidence intervals in 
equation (2.7) and (2.8), the two-sided confidence interval on r. 
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can be readily obtained. The 1-a two-sided confidence interval for 
is given by LJ <_ <. , where 


L * 


+ F o/2:n r - F 1 + (1 ' F a/2:n 1 ,- F a/2:n 1 ,n 3 )F a/2:n 1 ,n 3 F l 1 


rf 

n 1 + F a/2:n 1 r F a/2:n 1 ,n 2 F 2 


(2.9) 


-1 + F 


U X 


l-a/2:n,,- F l + (1 " F 1-l/2:n 1 ,- F l-a/2:n 1 ,r. 3 )F l-a/2:n 1 ,n 3 F l 1 


n 1 + F l-a/2:n r « F 1-o/2:n r n 2 F 2 


This will be proved in Chapter 3. 

If either limit is negative it Is replaced by iero. 

As mentioned in Chapter 1, due to the symmetry in the model. If 

a confidence interval fcr is obtained, the confidence interval 

for 6y can be easily obtained. The formulas for the confidence 

7 7 7 

Intervals for e T are given below where Oy = Ioy /(o^ +o £ )J. The 
1 -a lower confidence point for 6y Is 

-1 + F F + M-F"^ F IF f" 1 

a:n^,® 2 ' a:n 0 ,“ a:n 0 ,n, ; a:n-,n. 2 

L T (F,,F )-~-ri- 2 - -^^-(2.10) 

"2 ♦ F a:n 2 ,» F a:n 2 ,a, F l 


The 1-a upper confidence point for e T is 


* + F l-a:n 2 ,<» F 2 + ^’ F l-a:n ? ,» F l-a:n ? ,n 1 ^ F l-a:n ? ,n^ F 2 
u t (f,.f 2 )-*—- ■ -n --2-i- (2.D) 

* r i F. Fi 

2 l-a:n 2 ,~ l-a:n 2 ,n^ 1 






The 1-a two-sided confidence points for are 

-1 + F ,1 F, + (1 —F .1 F _ )F /9 _ F, 1 

a/ 2 :n 2 »® 2 a/Z.n^,® a/ 2 :n 2 »n 2 a/ 2 :n 2 ,02 2 

JJ - --T|-—- 

n ', + F ZO.„ F ,, F, 

and * a/ 2 .n 2 ,» 0 / 2 : 02 , 0 , 1 (? 

^ + F 1-a/2:n 2 »® F 2 ^ ^l-a/2:02,®^l-«/2:02,02^ 1-a/2:02,02 2 

L T " 7) 

n 2 + F 1 -a/ 2 :n 2 ,» F l- a / 2 :n 2 ,n 1 F l 

If any.limit is negative it is replaced with zero. 

From the confidence intervals obtained for and Oj » and by 
equation ( 2 . 1 ) we get the corresponding confidence intervals for p a 
and p t . 

The 1-a lower confidence interval for p a is 


0i ’ Ai rv 

where L A = L A (F,,F 2 ) which is defined in equation (2.7). 
The 1-a upper confidence intervals for p a is 


(J3* \ 


< P* < " 


where U A = u a^ f i ,F 2 ^ which is defined 1n equation (2.8). 
The 1-a two-sided confidence interval for p A is 

(*») = - • (A) 

where and are defined in equation (2.9). 

The 1-a lower confidence interval for py is 






17 


where Lj 
The 


where 

The 


Where U* 

( 2 . 12 ). 



= which Is defined In equation (2.10). 

1 -a upper confidence Interval for pj is 



Uy = Uy(F^,F 2 ) which Is defined In equation (2.11). 

1 -a two-sided confidence interval for o v is 
/ \ / \ 1 



and Ly are the confidence points given in the equation 
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3. Evaluation of the Methods 


In the preceding section confidence point(s) were determined 
for upper, lower and two-sided confidence intervals which have exact an 
confidence coefficient 1-a for specified restrictions, and for general 
conditions the confidence coefficient should be close to 1-a . In 

this section methods will be derived to determine how close 
the approximate confidence coefficients are to the specified ones. 

First consider the upper confidence interval 

PCL(F r F 2 ) < 0 A < -] = P L (e A ) (3.1) 

From equation (2.7) 


/ ,- F a:n, ,n > F „:n, .n F l’ 

I -!-' --!— i — if F, > F 


L(F 1 ,F 2 ) 


n, + F ' F F 9 

1 a:n^,® a:n^,n 2 2 


1 a. n ^, n^ 


If F, < F n 
1 a.n-|,n 3 


(3.2) 


From Section 1 
1 " 


F, ■ S, 2 /S, 2 for i * 1,2 ; also n^S^/e. « U i 


for i = 1,2,3 are Independent chi-square random variables with n. 
degrees of freedom respectively. 

Consider the following events 


A = {L(F r F 2 ) < e A ; F ] >0; F 2 > 0} 


■ < F 1 > F „: V ,y r 2 > "• 1 < F 1 .f 2 ) i 0) 
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B 2 * tF l ± F a: V n 3 ; F 2 > °* L < F l’ F 2 >i 0) * < 3 ' 3 > 

So B^n Bg * 0 and B^ U B fl where 

»■ Sfl.F 2 ,L(F,,F 2 )) :F,>0 ; F 2 >0iL(F r F 2 )>0} . 

Clearly P L (e A ) * P(A). We will write events B^ and B 2 in (3.3) as 

Bi * (F, > F n } ; B 9 = fF, < F „ } (3.4) 

1 1 a:n^,n 3 2 1 — a:n^,n 3 ' ' 


From this we obtain 

P L (8 A ) * P(A|B 1 )P(B 1 ) + P(A|B 2 )P(B 2 ) 

W ' PCL(F,.F 2 ) < e A |F, > F o;ni>n3 3F[F, > F a;n)i „ 3 ] 

+ P[L(F,.F 2 )<. A |F, iF a!v „HF, if. ;v „ 3 ] 


P L< 8 A>* P 


1+F 


:n.| ,"> F 1 + ^ F a:n.j ,® F a:n^ ,n 3 ^ F a:n.j ,n 3 F l 


-1 


-1 


n, + F "H F F, 

1 a:n^,« a:n^ ,n 2 2 


i e A 


F i >F „ 
I ol • n j 


* P I f l > ^n,.^ * i . a |F, < F a; „ )ir , 3 ]P[F, < F o: „ )t „ 3 ] 


But Is always greater or equal to zero so 
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and 


W ’ p 


niif v f 2 )l e A l Fi£f a: „ 1>n3 ]-1 


■ UF «I, ,- F a:n, ,n, )F (,:n, ,n/l' 


^—;-;-- 0 a 


n l + F a:n^,« F a:nj,n£ F 2 


F.>F „ 
1 o:n, ,n 




The first term in (3.5) gives 


P, ■ P 


- 1 +F <J ni > > F 1 *< 1 -F;.n 1 ,- F a;ni ,n 3 ) F a:n 1 ,n 3 F 1 
n, + F F " „ F 9 

1 o:n-j,<*> a:n-j,r2 

F a:n^ ,® +F l + ' F a:ni,» F a:n^ ,n 3 A:n-j ,n 3 F l 


< d. 


n.F + F . „ F, 

I a.,«> a:n-j ,n 3 2 


< 0 


p, = P 


•F S. 2 $, 2 *S, 4 *(F -F „ )F „ S, 
ain-j,® I o I a:n.j,°° a:n^,n 3 ' a:n^,n 3 3 


F,>F 

1 a:n^,n 3 


F I >F .n n 
I o.n-j.n^ 


n.F ~ $A 2 + F . n n S, 

1 q< n>|) 00 3 1 a* H| *^2 ^ • 


y - -6 - - < e. IF, >F <n n 

2^ 2 — AM a. n.j ,n 


Pi * p [-F $A. 2 + S, 4 + (F -F )F S 4 < 
1 L a:n-|,» 1 o 1 ' a:nj,® a:n^.n^ a:n^ ,n 3 3 — 


n.F S A A a + F , n n sAA.IF. > F . n 1 

1 a:nj,» 1 3 A Z 1 A'l a.n^,n 3 J 


p, .p 


WW t(F , )F /WW 

/n, 2 “'"I- °'V n 3 a "r"3 \ 6 3 2 / ' 
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(1-Y) 


-1 F a:n 1 ,n 3 + n 2 


n 2 p l 


| Y > a:n T n 3 


^ flop i 
a:n^ ,nj 2 1 


(3.8) 


Then (3.7) can be written as 
.P L (e A ) = [( C ! Y + C 2 0-Y) 2 Y _1 - C 3 (1-Y))x < U 2 ; 

Y ” F a:n 1 ,n/( F a:n 1 ,n 3 +n 2 p l ) ] 

+ P|Y < F n /(F „ +n oPl )l 

L — o:n|,n^ a:n-j,n 3 2 1 J 

= P [((C 1 +C 2 +C 3 )Y+C 2 Y~ 1 - (2C 2 +C 3 ))x < U 2 ; Y > a] + P[Y < A] 


where 


= P[g(Y)X < U 2 ; Y > A] + P[Y < A] 

9(Y) = (C^C^JY + C 2 Y _1 - (2C 2 +C 3 ) 

A = F „ /(F _ + n 0 p,) . 

a:n^,ny a:n^,n 3 2 1 


(3.9) 


With probability one Y e (0,1]. 

Consider the function 

(Cn+C/Cjy 2 - (2C ? +C )y + C ? 

g(y) - — 2 - - 3 —-—2-J-1 o < y < l 


g(y) = o iff (C 1 +C 2 +C 3 )y^ - (2C 2 +C 3 )y + C 2 = 0; 


Clearly 




25 


hence, g(y) has only two zeros for ye R, and they are 



(2C 2 +C 3 ) + [C 3 2 -4 Cl C 2 ] 1/2 

~~ 2lc7 + C 2+ c 3 l 


. <V C 3> - [ C 3 2 - 4C 1 C 2 3 1/2 

°2 " 2(C 1 +Cg+C 3 ) 


Also Urn g(y) - -® and Limg(y) = C, > 0. Since g(y) is continuous 
y-»0 + y-+l 

in (0,1] then at least one of the roots (6^,6 2 ) is in (0,1). But 


^jM = ( Cl +C 2 +C 3 ) - C^" 2 = C ] +C 3 +C 2 (l-y’ 2 ) , 0 < y < 1 

But for y e (0,1), ^1 - j< 0 . Also, (C^+C 3 ) > 0 and C 2 < 0 s 
hence, if ye (0,1) then 

> o . 
dy 

So g(y) is monotonic increasing in (0,1]; also since 

Urn g(y) = - ® and Lim g(y) = C, > 0 , only one of the two roots 
y-K) + y*l 

can be in (0,1]. This gives the following two possible forms for g(y). 



Figure 3.1 








Figure 3.2 


The number 6 is the root in (0,1); also Ae(0,l) so if g(x) < 0 
this implies that A < 6 . Next we show that g(x) < 0 for some value 


of A e (0,1). 


g(A) = (C 1 +C 2 +C 3 )a + C 2 A- ] - (2C 2 +C 3 ) 


( c l+ c 2 + c 3 )F a; 


^ F a:n,,n- +n 2 p l^ 


• (2C 2 +C 3 ) 


{ f a:n r n 3 + Vl) 2 F a:n r n 3 2 3 ' 

<Cl^2*C3>F| : n 1 .n/ C 2t F a:n,.n^ n 2 p l> 2 -< 2C 2^ C 3> F c:n 1 .n 3 < |r a :n,.n3* n 2 p l> 


a:n^»n^ ^ n 2 p ]* 


First note that the denominator of g(A) is always positive 
Denote the numerator by N so 

N - (C,*C 2 * C 3> F f:n,.n 3 * c J F f C 2 " 2 V + 2 C 2 Vl F ,:n, 

■ 2 C 2 F a:n 1 ,n 3 ' 2C 2 n 2 °l F ti:n r n 3 ' C 3 F «:n r n 3 ‘ C 3 n 2 t ’l F o:n, ,n 


. r 2 , r 2 ? r _ r 

C l F a:n,,n„ C 2 n 2 P 1 ‘ C 3VTa:«,,n, 
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Replace C^, C^. and by their values in (3.6) and obtain 


n 0 p iF 

2 1 q • n ^ * n ^ 


. n 2 2 C c „ 2 2 C 2 

Pi F i _ * Pi F _ 

2 1 a:n^ *» a:n-j ,n^ 2 1 arn^n^ 


ni P/j9* P _ « n iP ip a F _ r 

I 2 A ci«n| 1^2 3 1 2 A ctiDj )D 


n 9 ( @ A n i + ^) piF . n F 
7 A 1 1 A«i^i j® ot• ri| 

n*jP#}Q*F t 

3 2 A a.n^ 


L 2 p 2 F -n 2 P 2 F +F n ? 2 p 2 -nA^e.^+IJF If 

|^2 1 a:n, |» n 3 2 1 a:n ]» n 3 a*n«|»® 2 1 2 1 A 1 a:nj«n» 

HqPiPo^aF ,- n 
3 i 2 A a.n^ 


F . n F . n _ n 2 p, 2 (l-e.n 1 +l) 
a.n-j ,«> a.n^ ,n^ 2 1 A 1 

n 3 p l p 2 0 A F a:n lf n 2 


■ e fl n i n 9 2 Pl 2fr F _ 

A 1 2 1 a:n.j,® a:n-| »n^ 

n oP iP a F 
3 12 A a.n.| ,n 2 


< 0 


So the numerator is always negative which implies that g(x) will 
always be negative; and hence, A < 6 . 

Summarizing these results give 

a) g(y) i 0 iff y ^ 6 for y e (0,1) 

b) 0 < A < 5 < 1 ' ' (3.10) 

c) g(y) <0 iff y e (A,6) 

Using these facts to evaluate 

P[g(Y)X <. U 2 , Y > a] in (3.9) 





we obtain 


* P£g(Y)X < U 2 ; Y > A|A < Y < s] P[a < Y < 6] 

+ P[g(Y)X < U 2 ; Y > A|Y > 5] P[Y > 5] 

* 1 • P[A < Y < 5 ] + [g(Y)X < U £ ; Y > A; Y > 5 ] 

■ P[A < Y < 6] + P[g(Y)X < U 2 ; Y > 6] 

- Fy(6) - F y (A) + P [g(Y)X < U 2 ; Y > 5 ] 

where * 

Fy(y) = P[Y < y] 

From the original probability in (3.9) we get 

P L (0 A ) = F y ( 6 ) - Fy( a) + P[g(Y)X < U 2 ; Y > 6 ] + F y (a) 

“ F y ( 6) + P[g(Y)X < U 2 ; Y > 6] (3.11) 

Also 

P[g(Y)X < U 2 ; Y > 5] = P[g(Y)X + X < U 2 + X; Y > 6] 

■ p[g(Y) + 1 < U - 2 l X ; Y > s] 

c p[[g(Y) + i]" 1 ; Y > 6 ] (3.12) 

Let Z = X/(X+U 2 ) ; It follows that Z is distributed as Beta with 

/ r»i+ 1 I 3 n 2 . 

parameters \— 2 - ' ~Tl ' Also, ^ et n i be an even integer so that 

n.j/2 = for i = 1,2,3 and hence the are positive integers. 

If X.|, X 2 , X^ ’re independent chi-square random 
variables with , m 2 , degrees of freedom, respectively, and 
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'1 + X 2 


, . ,v^ 

2 a i + Xg + X 3 


then , Y ? are mutually Independent Beta random variables with 
parameters (m^/2, m^/Z) and ((m 1 +m 2 )/2,m 3 /2^ respectively. 

From the result in (3.12) we obtain (let v denote [g(y) +1]~*) 

p l ( 0 a ) - pftgm+ir 1 > Z; Y > «] + F y (6) 


Fy(6) + 


% y ^ 

/ 'r B (m, ,m_) m,-l m,-l m,+m,-l m..-l 

/ B Vizbh » ' «'*> ’ Z ' 3 <'-> 3 “ ^ 
6 0 J C 


1 «V -1 

r'r B (m,,m,) m,-l m~-1 m +m,-l 

■wy tovM? y °' y) z 

6 0 


m 2 -1 /m 9 -l\ rru-r-1 

Z v z r Jo) (-*) 2 ^y 

r=0 ' r ' 


m 1 1 rn.-r-l 

V B (m l’ m 3 )r(m 2 )(_1) 

F Y ( 6 ) + B (m 1 +m 3 ,m 2 ) r l r + V) r (m^rT 


Cg(y)+i] 


■1 


5 L 3 W/ 'J 

/y 1 " 1 '' 0-y)" 3 " 1 / z ( m l +m ? +m 3" r ‘^) d z dy 
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/ I _ 1 _ , -(m.+nu+nu-r-l) 

y 1 ( 1 -y) 3 j-taiilill ' 2 3 dy 


m2“r-1 


V 1 (_1) ‘ 

* Fy ^ + r5) "U J ^ 3 -1)! (m 2 -r-1;!r! (m 1 +n» 2 +m 3 -r- \) 


J m '’ 1 (l-y) 1 " 3 1 


* N/..\i.n I 


« [g(y) + U 




dy 


(3.13) 


From the expression In (3.13) w“ will compute p L ( 0 ^) 
Next consider the lower confidence Interval 

P[e A < u(f,.f 2 )] ■ Pu ( 6 A ) 

From equation (2.8) It follows that 

- FT 1 F, + (1-F -1 
l-a:ni ,•» 1 v 

U(F, 


(3.14) 


^ + F l-a:n],«* F l + ^ F l-a:nj ,» F l-a:n.| ,n^ F l-a:nj .n^l 

r r F 2 ) = ~~~n-7 ~ 

1 l-a:n-j l-a-n^ ,n 2 


if S, 2 /S, 2 > F, , n _ 
1 3 1- cufij, 


if S, 2 /S 2 < F, , n n 
1 3 I" cl n ^ 


U(F r F 2 ) * 0 

Substitute equation (3.15) into (3.14) to obtain 
Pu(e A ) * p[[g(Y) + l]' 1 <2; Y > «] 


(3.15) 


(3.16) 


where Y ~ Beta (m^.m^) , and Z ~ Beta (m^+m^,m 2 ) ; Y and Z 
are Independent. 
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Since g(y) is continuous in (0,1], then either both of the two roots 
of g(y) are in (0,1] , or neither of the roots are in the interval 
(0,1] , or the two roots are complex and not real. But A c (0,1] so 
that g(A) < 0 . For 

Fl * o: V n 3 


1-a:nj,n 3 + p 1 n 2 


(3.19) 


2 2 

_0 A P 1 n l n 2 F l-a:n 1 ,- F l-c.:n,,n, 

g( a) -- r - -— < o 

n 3 p 1 p 2 e A F l-a:n r n 2 


This follows from the proof of (b) in Equation (3.10). 

From the above it Is clear that the two roots of g(y) are in 
(0,1], which results in the following form for g(y). 



From equation (3.16) we evaluate Py(e^) as Follows 

_ , nu-r-l 

V' (-1) 2 (m 1 +m 2 +m 3 -l)! 

P i/ e A^ * ^ TmpTJ! (m^-l)! (mg-r-l)TrV (m^^^+m^r-Ty 


i 
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► 

7 1 hIt -1 m ,-1 

y 1 0-y) 3 

« 

The two-sided confidence Interval Is 

P[L 1 (F r F 2 ) < e A < U 1 (F 1 ,F 2 )3 « P 2 (e A ) (3.21) 

where Is given by (3.2) with o replaced by a/2 and 

U] (F].F^) given by (3.15) with 1-a replaced by l-a/2 . 

To show (3.21) let A and B be the following events 

A - {L 1 (F 1 ,F 2 ) < 6 a ) B - {e A < U^F^Fg)) 

Then 

P 2 (« A ) * P(A f\ B) (3.22) 

From (3.22) we obtain 

' p < A fi B ) * 1 - p(rrii) -1 - p(aUb) 

: 1 ~ [P(K) + P(B) - P(AflB)] 

■ 1 - P(ff) - P(ff) + p(ffOff) (3.23) 


[g(y)+i] 


(mj+m2 +m 3‘ r "U 


- + 1 


dy 


Also 


P(^) « 1 - P L (e A ). P(B) - 1 - P u (e A ) 
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P(*HB) « p(e A < L(F r F 2 ); U(F r F 2 ) <D A )=0 

\ 

If L(F^,F 2 ) iU(Fj,F 2 ) with probability one. Next we will show 
that this is indeed the case for n-j > ? , n 2 >_ 2 and a <_ .10 . Vie 


must prove 


-1 + F*)_ r, + (1-F”] F ) F n f: 1 

aini,® i a:n.| *® a:n^ .n^' oin^.n^ 1 

r.i ♦ F-L FT „ F, 

1 a:n«| a:n^ • 2 

-1 + ^ + (l-p; 1 Pi n )Pl n n 

l-a:n^i«- 1 l-a:n^,*» l-arn^n^ l-ain-j.n^ 1 

n l + F l-o:r. r « F l-a:n r n 2 F 2 


which we write as 


•1 + A-jFj + BjFj 1 -1 + AgFj + B^ -1 


’1 H r 2 


n l + C 2 F 2 


This implies 


U(F r F 2 ) < L(F r F 2 ) 


To prove equation (3.24) we use Lemma l, 2, and 3. (See Ghosh (1973)). 


Lenina l . F_ „ Is decreasing in n 0 if n, > 2 , n 0 > 2, 

- o.n 1 ,n 2 c i c. — 


a < .10 


Lemma 2 . F _ is decreasing in n, if n. > 2, n, > 2, 
.. i... _ a.n^ ,n 2 i i c 


a < .10. 
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Lenina 3 . F i_ a;n , n Is Increasing in n 2 if ^ > 2 , 
n 2 > 2, a < .10 . 

Proof: 

F a:nf,n 2 < F a:n r n 2 1f n 1 > n l to Len ™ 2 > • 


So 


F' 1 < F’ 1 

l-a:n 2 ,n^ l~on 2 »n^ 


And from this the result follows. 

Using the Lemmas in expression (3.24) we have 

1) C 2 < C<| since 


then 


l-a:n^ ,«• l-a:n^,n 2 — F o:n^ i® F o:n^ ,n 2 


-1 + a^, + b^f: 1 . i + a,f, + b,f: 

r l r ^ i *1 • I 11 11 

C 1 ~ 1 ’ C 2 - 1 S ° n, + 


-1 


2) F"] < Fi] 50 A 1 < 

a.n-j » I I — 


and 


-1 + A ] F 1 + B^ 1 - 1 + A„F, + BiF 


-1 


n l + C 1 F 2 


2 i °ri 


n 2 + C 2 F 2 


so 


3) F' 1 F > F”^ F 

otn^*** a:n^ .n^ — l-a:n^»* l-a:nj»n.. 


F a:nj .^arn^ ,n^ - F l-a:n.j ,® F l-a:n^» 


(1-F F ) < (1-F: 1 F. ) 

a:n.j,« a:n^ ,n^ - l-a:n^,« l-«:n 1 ,n 3 ' 


and 
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--1 


»«» F a:nj »n,^ F a:nj»n^ — ^ F l-a:nj ,®^l-a:n'j ,n^ l-am^.n^ V 


then 

and 


B 1 — B 2 


- 1 + A,F, + B^f; 1 - 1 + A 2 F 1 + B/’ 1 


Yi °Y 1 
"l + C 1 F 2 


n l + C 2 F 2 


Going back to (3.23) we get 

P 2 (b a ) * 1 - ('-P L (»A>) • (’-w) * 0 

P 2 (9 A ! - P L (e A ) * P U (« A ) - 1 (3.25) 


Numerical Integration will be used to evaluate the integrals 
In P L 0 A ) , P^e^), and P2^ 0 a^’ ^ he method usecJ sim P son ' s Rule 
which gives 

f\y)dy £fu) + 4 ^ a+ h) + 2f(a+2h) 

+ 4f(a+3h) + 2f(a+4h) + ... + 2f(a+(n-2)h) + 4f^a+(n-l)hj 

+ f(b)] 

where h Is the length of the subintervals In which (a»b) has been 
partitioned. We used h = .001. 

Various combinations of sample sizes were cnosen, and for each 
of the combinations several sets of , and , and different 
a's were used. 
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A list of the values of I, J, p^, p^, and a used Is given 

below 

(I.J) • (3»3);(3,7);(3,11);(5,3);(5,5);(5,11);(7,3); 

(7»7);(7,11);(9,3);(9,9);(9,11);(11,3);(11,7);(11,11) 
o a ? /o £ 2 : .003, 1, 10, 100 (4.1) 

o 2 /c 2 : 0, 1, 10, 100 
l e 

1-a : .90, .95, .975, .99 

1 + J °k ,v t 
p2 • 1 + I 0 T^/° e ^ 

The results are given in Tables 2, 3, and 4. The first five columns 
in each table show the I, J, n^, and n^ values. Each of the remaining 
columns give the results for the specified 1-a confidence coefficients. 
For each of these columns the ranges that the confidence coefficients vary 
(as p.j and p£ vary over their allowable values) are shown. 






The approximate confidence coefficients corresponding to the 
lower confidence intervals are very close to the specified ones, shown 
in Table 2. The approximate confidence coefficients for the upper 
confidence intervals range from values very close to the specified ones 
to some high values, shown in Table 3. Similar results hold for the 
two-sided confidence limits shown in Table 4. 

For I fixed the range of the confidence coefficients is closer to 
the nominal value when J increases. For J fixed the upper bound of 
the range is higher when I increases. So the confidence coefficient 
will be better when J is large. However, in general the approximate 
confidence procedures seem adequate if I and J are each 5 or bigger. 






Ranges of confidence coefficients for lower confidence Intervals on 
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Table 4 

Range of confidence coefficients for a two-sided confidence 
interval on 


I 

J 

n l 

n 2 

n 3 

1 - 

a 

= .90 

1 - 

a 

O) 

ii 

3 

3 

2 

O 

4 

.90103 

- 

.93208 

.95072 

- _ 

.97000 

3 

7 

2 

6 

12 

.90043 

- 

.91349 

.95030 

- 

.96000 

3 

11 

2 

10 

20 

.89950 

- 

.90849 

.94827 

- 

.95702 

5 

3 

4 

2 

8 

.90096 

- 

.93896 

.95062 

- 

.97273 

5 % 

5 

4 

4 

16 

.90077 

- 

.92684 

.95051 

- 

.96725 

5 

11 

4 

10 

40 

.89929 

- 

.91272 

.95031 

- 

.95884 

7 

3 

6 

2 

12 

.90094 

- 

.94310 

.95604 

- 

.97436 

7 

7 

6 

6 

36 

.89871 

- 

.92475 

.94984 

- 

.96607 

7 

11 

6 

10 

60 

.89896 

- 

.91673 

.94955 

- 

.96130 

9 

3 

8 

2 

16 

.90092 

- 

.94619 

.95059 

- 

.97647 

9 

9 

8 

8 

64 

.89835 

- 

.92374 

.94929 

- 

.96546 

9 

11 

8 

10 

80 

.89878 

- 

.92016 

.94938 

- 

.96334 

11 

3 

10 

2 

20 

.90090 

- 

.94863 

.95058 

- 

.97782 

11 

7 

10 

6 

60 

.89790 

- 

.93203 

.94926 

- 

.97006 

11 

11 

10 

10 

100 

.89868 

- 

.92315 

.94927 

- 

.96510 
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